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I.  INTRODUCTION 


Because  of  their  unique  properties,  brittle  materials  must  be  used  in 
many  military  applications.  Such  materials  have  superior  strength  at  high 
temperatures  and  are  therefore  used  in  reentry  nosetips,  e.g. , graphite  and 
beryllium.  Investigation  of  the  feasibility  of  using  silicon  nitride  and 
silicon  carbide  for  turbine  motors  is  under  way  by  DARPA.  Also,  virtually 
all  materials  that  are  transparent  in  either  visible  or  infrared  wavelengths 
are  brittle.  These,  therefore,  must  be  used  in  radomes  and  laser  windows, 
e.g.,  quartz,  silicon,  zinc  selenide,  and  cadmium  telluride. 

Brittle  materials  usually  fail  by  fracture.  Unlike  ductile  materials, 
they  exhibit  a wide  dispersion  in  fracture  stress,  and  this  dispersion  must 
be  taken  into  account  in  design.  This  requires  the  use  of  fracture  statistics, 
which  are  currently  inadequately  understood. 

Fracture  statistics  can  be  obtained  in  the  laboratory  for  simple  tension 
and  for  some  biaxial  loading  conditions.  However,  the  stresses  encountered 
in  service  are  generally  triaxial.  These  cannot  be  reproduced  in  the  laboratory, 
and  even  if  they  could,  the  need  for  statistical  data  at  each  of  a large  number 
of  combinations  of  stress  w__ld  make  testing  prohibitively  expensive.  Thus, 
there  is  a need  for  a theory  that  shows  how  to  use  laboratory  data  in  simple 
tension  to  predict  the  fracture  statistics  associated  with  any  combination  of 
stresses. 

Until  very  recently,  the  best  tool  to  use  was  Weibull's  weakest  link 
theory  (R.ef.  1).  This  is  a statistical  theory  that  assumes  that  fracture  is 
caused  by  the  presence  of  flaws  of  unspecified  nature.  In  this  theory,  the 
fracture  data  in  simple  tension  are  represented  approximately  by  proper 
choice  of  three  parameters  in  a formula  given  by  Weibull.  For  some  cases, 
a reasonable  fit  can  be  obtained  with  a two -parameter  formula  obtained  by 
making  one  of  the  parameters,  (the  minimum  stress  that  can  cause  frac- 
ture), equal  to  zero.  Weibull  gave  a rule  for  deducing  polyaxial  stress 
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statistics  from  uniaxial  and  showed  how  to  apply  it  in  the  simpler  two- 
parameter  formula.  However,  he  did  not  treat  the  usual  case  where  ou  =£  0, 
as  in  the  case  of  graphite.  Moreover,  the  rule  for  converting  from  uniaxial  to 
polyaxial  statistics  must  be  regarded  as  intuitively  plausible  rather  than 
theoretically  or  experimentally  demonstrated.  A very  approximate  technique 
for  extending  Weibull's  approach  to  the  case  au  ^ 0 was  proposed  by  Dukes 
(Ref.  2),  but  no  experimental  check  on  its  accuracy  was  attempted.  A major 
shortcoming  in  Dukes'  approach  is  that  his  treatment  is  internally  inconsistent. 

In  1970,  a reexamination  of  the  fundamentals  of  fracture  statistics  was 
undertaken  by  Batdorf.  His  approach  was  a weakest-link  theory  in  which 
extreme  value  statistics  were  applied,  but  it  differed  from  Weibull's  in  two 
major  respects.  First,  the  flaws  were  assumed  to  be  cracks,  with  the 
directional  sensitivity  of  cracks  to  the  applied  stresses.  Second,  rather  than 
a three -parameter  approximate  fit  to  the  data,  a Taylor  series  expansion 
was  used  to  fit  the  data  to  arbitrary  precision.  The  theory  was  applied  to 
uniaxial  and  biaxial  test  data  on  POCO  graphite  and  found  to  be  in  very  good 
agreement  (Ref.  3). 

Unlike  POCO,  the  advanced  graphites  of  greatest  interest  for  nosetip 
application  are  not  isotropic,  but  exhibit  a moderate  degree  of  anisotropy. 
Therefore,  the  theory  was  extended  to  cover  this  case  (Ref.  4).  The  theory 
has  also  been  extended  to  cover  the  case  of  materials  with  surface -distributed 
cracks  (Ref.  5).  Application  of  the  theory  to  extensive  data  on  pyrex  cylinders 
indicates  that  it  is  in  better  agreement  with  experiment  than  any  other  available 
theory. 

The  preceding  theory  permits  the  use  of  fracture  data  in  simple  tension 
in  order  to  determine  the  probability  of  failure  of  any  volume  of  material  in 
an  arbitrary  stress  state.  It  can  therefore  be  incorporated  into  a finite  - 
element  computer  code  for  determination  of  the  probability  of  survival  of  a 
nosetip  during  reentry. 

An  ideal  statistical  theory  of  fracture  of  brittle  materials  must  take 
three  things  properly  into  account:  (1)  extreme  value  statistics,  (2)  fracture 
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mechanics,  and  (3)  material  properties  and  micro  structure.  Weibull's  theory 
(Ref.  i)  satisfies  the  first  requirement  approximately,  but  neither  of  the  others. 
The  later  work  described  previously  (Refs.  3 through  5)  incorporates  the  first 
two  requirements;  but  not  the  third.  Very  recently,  efforts  have  been.under- 
taken  to  develop  statistical  theories  directed  toward  satisfying  all  three 
requirements. 

Theories  that  do  not  satisfy  the  third  requirement  have  to  depend  on 
experimental  data  in  one  stress  state,  e.  g. , simple  tension,  ‘.n  order  to 
establish  the  distribution  of  crack  strengths;  fracture  statistics  can  then  be 
determined  for  other  stress  states.  The  objective  of  theories  that  satisfy  all 
three  requirements  is  to  predict  the  distribution  of  crack  sizes  from  material 
microstructures.  Such  theories  have  many  advantages:  they  require  much 
less  experimental  data,  they  can  be  used  to  determine  the  changes  in  material 
processing  needed  for  desired  improvements  in  fracture  behavior,  and  they 
should  permit  reliable  extrapolation  of  the  curve  of  failure  probability  versus 
stress  to  very  low  probabilities  of  failure. 

Progress  in  the  development  of  theories  that  satisfy  all  three  require- 
ments is  being  made  in  two  areas.  One  is  the  statistical  theory  of  fracture  for 
brittle  materials  with  preexisting  integranular  cracks.  McClintoclc  (Ref.  6) 
developed  a two-dimensional  model  for  such  a material  in  which  line  cracks 
were  all  parallel  and  normal  to  the  direction  of  tension.  The  corresponding 
three-dimensional  theory  for  uniformly  distributed  penny-shaped  cracks  was 
developed  by  Batdorf  (Ref.  7).  This  theory  should  apply  to  any  polycrystalline 
ceramic.  The  second  area  is  the  theory  of  the  initiation  and  growth  of  cracks 
in  materials  such  as  graphite.  Such  a theory,  which  predicts  the  stress -strain 
relation  as  well  as  the  fracture  stress,  has  been  proposed  by  Buch,  Zimmer, 
and  Meyer  (Ref.  8).  In  both  cases,  the  theories  should  be  extended  and  refined 
in  order  to  improve  their  accuracy  and  to  make  them  applicable  to  polyaxial 
stress  conditions. 

The  prime  objective  of  the  present  work  is  to  further  develop,  refine, 
compare  with  experiment,  and  modify  the  theories  as  necessary  - those  in  which 
the  critical  stress  frequency  distribution  of  cracks  is  determined  experimentally. 
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and  those  in  which  the  distribution  is  based  on  theoretical  considerations.  For 
convenience,  the  two  classes  of  theory  are  referred  to  herein  as  experimentally 
based  and  theoretically  based,  respectively.  The  experimentally  based 
theories  will  be  of  greater  use  in  structural  design  in  the  short  run.  The 
theoretically  based  approaches  have  more  ultimate  potential  and  are  of 
greater  value  to  materials  producers.  They  are  also  better  when  large 
changes  in  volume  or  failure  probability  are  required  in  structural  design. 

In  the  proposal,  the  tasks  that  were  specifically  identified  were  the 
following: 

1.  . Experimentally  Based  Approach 

a.  Investigate  the  effects  of  porosity  on  fracture  statistics. 

If  possible,  include  porosity  as  a parameter  in  the  theory. 

b.  Reformulate  the  theory  to  include  the  statistical  conse- 
quences of  material  variability. 

c.  Extend  the  theory  to  include  also  predominantly  com- 
pressive states  of  stress. 

d.  Conduct  an  extensive  literature  search  for  fracture 
data  with  special  emphasis  on  graphite.  Check  the 
validity  of  experiment  by  comparison  with  these  data. 

Modify  the  theory  if  necessary. 

2.  Theoretically  Based  Approach 

a.  Refine  the  theories  for  greater  accuracy  and  to  permit 
use  of  more  realistic  material  models,  e.g. , variable 
grain  size,  porosity,  and  directional  orientation. 

b.  Extend  the  theories  to  apply  to  polyaxial  stress  conditions. 

c.  Compare  the  theories  with  experiment  and  revise  as 
necessary. 


II.  FY  75  ACCOMPLISHMENTS 


A.  LITERATURE  SEARCH 


There  are  two  reasons  why  a literature  search  for  fracture  data  is 
significant  to  this  investigation: 

1.  By  artful  selection  of  data,  or  even  innocent  accidental 
choice,  almost  any  theory  can  be  corroborated  or 
discredited.  Only  after  an  extensive  survey  can  the 
experimental  facts  be  stated  with  confidence. 

2.  The  theory  is  fairly  complicated,  and  refinements  will 
probably  be  more  complicated.  It  is  not  desirable  to 

put  too  much  effort  into  refinements  until  it  is  determined 
if  the  basic  treatment  is  in  accord  with  data. 

The  literature  search  for  data  and  analyses  of  fracture  statistics  of 
brittle  materials  is  progressing  satisfactory.  Western  Research  Contract 
Center  has  completed  its  computerized  search  for  appropriate  references. 
This  search  covered  the  NASA  Data  Bank  Abstracts  from  19  69  to  the  present. 
Chemical  Abstracts  from  1970,  and  Engineering  Abstracts  back  to  1971. 
Earlier  sources  are  not  in  the  computer  banks  and  these  are  being  located 
in  the  time -honored  manual  fashion  by  M.  Adams  at  the  Univ,  of  California, 
Los  Angeles.  He  is  digesting  all  relevant  papers  and  reproducing  appro- 
priate graphs  and  tables  for  ready  reference.  This  work  is  under  way  and 
will  be  completed  before  the  end  of  FY  75.  A sample  page  of  the  computer 
search  report  and  a sample  digest  of  a paper  are  presented  in  Appendix  A. 

B.  FRACTURE  STATISTICS  OF  BRITTLE  MATERIALS 
WITH  SURFACE  CRACKS  (Ref.  5) 


The  most  important  differences  between  this  theory  and  the  earlier 
one  with  volume -distributed  cracks  (Ref.  3)  lie  in  the  assumption  that  the 
surface  crack  planes  are  normal  to  the  surface.  Thus,  the  probability  of 
failures  in  simple  tension  takes  quite  different  forms  for  the  two  cases 
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(1) 


Pf(a)  = 1 


Pf(a)  = i 


dN(a  J . . ] 

"d5 I1  “ Jacv/(r  Kr. 

cr 

dN(a  ) . 

— ; — — cos"  la /or  da 
dorcr  V cr  cr. 


(2) 


where  Eqs.  (i)  and  (2)  apply  to  volume  and  surface -distributed  cracks, 
respectively.  In  both  theories,  the  function  N(acr),  which  represents  the 
number  of  cracks  per  unit  volume  (area)  with  a critical  stress  equal  to  or 
less  than  a^,  is  found  by  matching  uniaxial  test  data.  Knowledge  of  dN/dacr 
permits  evaluation  of  corresponding  expressions  for  probability  of  fracture 
under  an  arbitrary  polyaxial  state  £ by  using 


Pf(cr)  = 1 


P £(o)  = 1 


(3) 


(4) 


where  0 is  a solid  angle  that  corresponds  to  orientations  in  space,  and  u>  is  a 
planar  angle  that  corresponds  to  orientations  in  a plane. 

In  Ref.  3,  the  technique  used  to  solve  Eq.  (i)  for  dN/d a was  to  expand 

v»  A 

it  as  a Taylor  series  and  evaluate  the  coefficients  so  as  to  fit  the  fracture  data 
in  simple  tension  at  a selected  number  of  stress  levels.  This  required  trun- 
cating the  series  and  inverting  a matrix.  It  has  subsequently  been  discovered 

that  dN/da  can  be  obtained  from  either  Eq.  (I)  or  Eq.  (2)  by  solving  them  as 
c r 

integral  equations.  This  avoids  truncation  errors  as  well  as  some  danger  of 
numerical  instabilities  and  is  therefore  preferable.  The  solution  of  Eq.  (i) 
is  given  in  Appendix  B. 
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The  corresponding  solution  of  Eq.  (2)  has  also  been  developed  and 

incorporated  into  a revised  version  of  Re.f.  5,  which  will  be  published  later. 

However,  a very  surprising  fact  was  recently  discovered:  for  given  uniaxial 

failure  statistics  Pf(a  , o'),  the  predicted  biaxial  fracture  statistics  P,(cr  ,o  ) 
i x x x y 

are  the  same  when  either  Eqs.  (1)  and  (3)  or  Eqs.  (2)  and  (4)  are  used;  i.  e. , 
the  volum-  theory  and  surface  theory  lead  to  the  same  polyaxial  stress  failure 
result,  although,  of  course,  the  functions  VN(a  ) and  AN(acr)  are  very  dif- 
ferent. Thus,  for  the  sake  of  economy  of  analytical  techniques  and  computer 
programs,  the  report  is  being  revised  to  exploit  the  already-developed  volume 
theory  approach.  This  will  not  change  the  numerical  results,  but  will  con- 
stitute quite  a change  in  the  point  of  view  and  the  depth  of  understanding  of 
the  problem. 

C.  GRAPHITE  MODEL 

The  theories  discussed  in  Refs.  1 through  7 assume  a brittle  material 
with  preexisting  randomly  oriented  cracks  that  remain  unchanged  during  the 
loading  process  until  the  critical  stress  of  the  weakest  crack  is  exceeded. 

At  that  time,  the  crack  becomes  unstable  and  grows  without  limit;  as  a result, 
the  specimen  fails. 

An  ingenious  material  model  for  graphite  behavior  has  been  proposed 
(Refs.  8 and  9)  in  which  the  cracks  are  visualized  as  being  created  by  the 
loading  process.  This  theory  considers  each  grain  in  a typical  commercial 
graphite  as  a crystal  and  exploits  the  fact  that  a crystal  of  graphite  is  very 
weak  in  the  c-direction.  Thus,  random  aggregations  of  grains  that  all  have 
approximately  the  same  c-orientation  will  easily  fracture  when  a tensile  stress 
is  applied  in  that  direction.  An  increase  in  stress  will  result  in  a relaxation 
in  the  degree  of  orientation  needed  for  crack -like  behavior  and,  therefore, 
more  and  larger  cracks.  The  nonlinearity  in  the  stress -strain  curve  is 
attributed  to  the  increasing  compliance  as  the  cracks  increase  in  number  and 
size.  Fracture  of  the  material  occurs  when  the  radius  r of  the  largest  crack 
exceeds  «.he  critical  size  for  the  applied  stress  as  given  by  the  fracture 
mechanics  equation 
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The  general  concepts  as  outlined  in  the  preceding  paragraph  are 
probably  basically  correct.  Even  though  conceptually  simple,  a precise 
treatment  of  this  model  is  mathematically  intractable,  and  approximations 
are  required.  The  theory  has  been  successful  in  accounting  approximately 
for  the  uniaxial  stress-strain  and  fracture  characteristics  of  typical  graphites. 
In  an  effort  to  improve  the  accuracy,  some  work  has  been  done  in  the  present 
investigation.  The  main  innovations  are  (1)  the  derivation  of  analytical  exr 
pressions  for  results,  where  practicable,  instead  of  depending  on  numerical 
analyses  by  computers,  and  (2)  a revised  criterion  for  crack  size.  In  the 
theory  of  Refs.  8 and  9,  an  oriented  penny-shaped  array  of  N grains  behaves 
like  a crack  only  when  the  component  of  the  applied  stress  in  the  c direction 
of  every  grain  in  the  array  exceeds  the  critical  stress  of  the  grain.  In  effect, 
this  gives  an  array  a strength  equal  to  that  of  the  most  misaligned  grain.  In 
the  present  investigation,  it  is  assumed  that  the  strength  of  the  array  is  equal 
to  the  average  strength  of  the  constituent  grains.  An  analysis  based  on  this 
revised  assumption  is  given  in  Appendix  C.  This  analysis  includes  a treat- 
ment of  porosity  for  the  case  where  pores  are  small,  i.  e. , approximately 
grain-sized. 

The  revision  has  resulted  in  improved  agreement  with  experiment,  but 
there  are  some  uncertainties  in  the  legitimacy  of  some  of  the  approximations 
used  that  need  to  be  cleared  up  before  publication.  One  concerns  the  elastic 
modulus  of  a material  with  uniformly  distributed,  partially  oriented  cracks 
such  as  those  implied  by  the  theory.  Another  is  the  legitimacy  of  the  neglect 
of  crack  interaction,  without  which  the  analysis  becomes  exceedingly  complex. 
The  latter  question  is  discussed  in  Section  IVE. 

D.  APPROXIMATE  METHODS  OF  DETERMINING  THE  FRACTURE 
STATISTICS  FOR  POLYAXIAL  LOADING 

In  his  original  treatment  of  the  statistics  of  fracture,  Weibull  (Ref.  1) 
introduced  a two-parameter  functional  form  for  the  relationship  between 
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, simple  tension  and  probability  of  fracture.  He  also  showed  how  to  compute 

the  corresponding  statistics  of  failure  for  bending,  torsion,  and  other  stress 
states  that  involve  only  uniaxial  tensile  stresses.  These  applications  have 
become  well  known  and  are  widely  used. 

In  problems  that  involve  biaxial  or  triaxial  tensile  stresses,  the  situa- 
tion is  more  complicated.  Weibull  gave,  without  formal  proof,  a procedure 
for  treating  such  problems  and  showed  how  to  apply  it  in  some  simple  cases. 

Some  investigators  have  expressed  doubts  concerning  the  rigor  of  Weibull's 
treatment  of  polyaxial  stress  states,  arid  there  are  indications  that  these 
doubts  were  later  shared  by  Weibull  himself.  Since,  in  addition,  rather 
tedious  calculations  are  required  for  each  polyaxial  stress  state,  in  practical 
structures  that  involve  continuously  varying  stress  states,  there  is  a natural 
tendency  to  use  approximations.  One  simple  approximation,  which  constitutes, 
in  fact,  the  only  technique  suggested  for  handling  polyaxial  stress  statistics 
" in  a recent  treatise  on  fracture  (Ref.  10)  is  to  assume  that 

| ■ WPs(a2>PS(o3>  (6) 

i 

| 

5 

where  Pg  is  the  probability  of  survival,  and  ci»<r2'tr3  are  principal  stresses. 

j The  analysis  in  Appendix  D points  out  that  the  above  approximation  is  often 

; seriously  unconservative.  An  exact*  treatment  is  also  given  of  failure  statis- 

| tics  in  biaxial  tension  for  cases  where  the  uniaxial  tensile  properties  obey 

j Weibull's  two -parameter  formula.  In  addition,  a conservative  approximation 

is  given  that  is  quite  accurate  in  many  cases  of  practical  interest. 

One  application  of  the  conservative  approximation  is  the  failure  statistics  | 
I of  laterally  loaded  Zn  Se  disks.  The  data  in  question  were  obtained  by  Univ. 

of  Dayton  Research  Institute  (UDRI)  as  part  of  the  Air -Borne  Laser  Laboratory 
‘ • program  at  AFWL.  In  Fig.  1,  the  unconservative  nature  of  the  prediction 


■ * 1 

Exact  from  the  standpoint  of  theories  in  Refs.  1 and  3. 
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RELIABILITY, 


based  on  Eq,  (6)  is  confirmed.  The  approximation  proposed  in  Appendix  D 
appears  to  be  much  better. 

E.  .STATISTICAL  MODELS  FOR  BRITTLE  FRACTURE  (R.  H.  Huddlestone 

and  D.  C.  Pridmore- Brown) 

*■-  "] 

Two  aspects  of  the  statistics  of  brittle  fracture  have  been  studied:  (1) 
the  relationship  between  extreme -value  statistics  and  possible  fracture  models, 
and  (2)  the  combinatorial  problem  of  calculating  the  flaw -size  distribution  in 
a material  specimen  that  is  assumed  to  consist  of  a regular  array  pf  identical 
grains,  each  with  an  equal  and  independent  probability  of  being  cracked.  The 
status  of  this  work  is  summarized  below.  Derivation  of  the  crack-size 
distribution  is  discussed  in  more  detail  in  Appendix  E. 

According  to  Gumbel  (Ref.  li),  the  first  application  of  the  asymptotic 
theory  of  extreme  values  to  fracture  problems  was  made  by  Epstein  (Ref.  12). 
In  this  and  later  work,  the  assumption  was  made  that  the  breaking  strength  of 
the  material  sample  is  a linear  function  of  the  size  of  the  largest  flaw  alone. 
This  assumption,  which  is  claimed  to  have  gone  unnoticed  by  other  workers, 
according  to  McClintock,  was  justifiably  criticized  by  him  in  a paper  (Ref.  6) 
in  which  he  developed  an  idealized  model  that  does  not  correspond  to  an 
asymptotic  distribution  of  extreme -value  theory.  McClintock1  s position  was 
further  corroborated  by  Batdorf  (Ref.  7)  with  a more  realistic  three- 
dimensional  model.  Here,  we  try  to  clarify  this  situation  and  analyze  the 
relevance  of  extreme -value  statistics  in  this  context. 

If  it  is  assumed  that  a material  specimen  under  tension  fractures  at  its 
weakest  cross  section,  then  a fracture  test  corresponds  directly  to  selecting 
that  cross  section  of  smallest  breaking  strength  from  the  large  sample  of 
cross  sections  in  the  specimen.  In  other  words,  the  breaking  strength  of  the 
specimen  can  be  expected  to  be  an  extreme-value  statistic  irrespective  of  the 
inadequacy  of  the  above  linearity  assumption. 

If  the  breaking  strength  exhibits  an  asymptotically  stable  distribution 
as  the  length  of  the  test  specimens  is  increased,  then,  necessarily,  by 
Gn'cdenko's  theorem  of  order  statistics,  the  distribution  will  be  the  third 
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asymptotic  one.  It  should  be  realized  that  this  distribution  was  introduced  on 
a purely  phenomenological  basis  by  Weibull  (Ref.  1)  and,  consequently,  is 
frequently  referred  to  as  the  Weibull. distribution.  However,  it  should  be 
emphasized  that  the  present  conclusidxgfa'and  their  justification  are  not  based 
on  Weibull's  model,  but  rather  on  a theorem  of  mathematical  statistics. 


There  appears  to  be  no  compelling  reason  why  the  distribution  of 
breaking  strengths  should  be  asymptotically  stable  on  physical  grounds.  It 
is  this  point  that  dictates  the  applicability  of  the  statistical  conclusions  stated 
above.  Gumbel  and  Gn'edenko  in  fact  give  examples  of  extreme-order 
statistics  that  are  not  asymptotically  stable  and  whose  distribution  does  not 
approach  any  of  the  three  asymptotic  laws.  The  models  of  McClintock  and 
of  Batdorf  are  both  of  this  type,  and  the  rescaled  shape  of  the  distribution 
changes  with  the  size  of  the  specimen. 

F.  FLAW -SIZE  MODEL  BASED  ON  THEORY  OF  RUNS  (R.  H.  Huddlestone 
and  D.  C.  Pridm  ore -Brown) 


The  work  on  the  flaw-size  distribution  is  an  effort  to  develop  a mathe- 
matical model  that  will  permit  an  explicit  treatment  and,  therefore,  a quanti- 
tative assessment  of  the  errors  implicit  in  the  more  heuristic  models  that 
are  currently  used.  Preliminary  theoretical  analysis,  supported  by  small 
specimen  calculations  made  with  the  Aerospace  on-line  APL  system,  con- 
vinced us  that  many  of  the  approximate  estimates  of  flaw-sir-e  distribution 
in  current  use  may  involve  large  errors  in  certain  applications. 

We  eventually  noticed  that  the  theory  of  runs,  a well -developed  branch 
of  probability  theory,  would  permit  a direct  calculation  for  the  highly  over- 
idealized  one -dimensional  case.  Such  a prototype  could  then  be  compared 
with  the  heuristic  formulas  at  least  in  one  dimension  and  could  prossibly  be 
used  to  improve  two-  and  three-dimensional  intuition  as  well.  In  any  weakest 
link  theory  that  relates  strength  to  crack  size,  statistical  questions  concerning 
the  distribution  of  crack  size  become  important.  Such  questions  are  difficult 
to  answer  in  general,  and  it  has  been  common  practice  to  make  various 
simplifying  assumptions.  One  is  the  assumption  that  crack  interaction  is 
negligible.  Intiutively,  this  assumption  appears  valid  in  the  limit  where  the 
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probability  of  finding  a crack  of  a given  size  is  very  small,  but  it  is  not  clear 
just  when  it  breaks  down. 
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In  order  to  study  this  question,  a simple  model  is  considered  in  order 
that  it  can  be  treated  analytically  without  neglecting  crack  interactions.  In 
this  model,  the  material  is  represented  by  a single  row  of  grains,  each  with 
a given  probability  p of  being  cracked.  A crack  is  defined  as  a row  of  one  or 
more  contiguous  cracked  grains,  and  we  consider  the  probability  that,  in  a 
sample  of  size  n,  i.e.,  in  a row  of  n grains,  there  exists  one  or  more  cracks 
of  size  J or  greater. 

The  mathematical  theory  that  permits  conclusions  merely  from  the 
order  in  which  the  elements  of  the  sample  appear  is  called  the  theory  of  runs 
and  is  treated  in  books  on  statistics.  In  particular,  some  progress  toward 
answering  the  above  question  was  made  by  Mood  (Ref.  13).  However,  he 
calculated  the  joint  probability  distribution  for  all  crack  sizes  less  than  J 
and  the  total  number  of  cracks  2J,  Thus,  for  even  moderate  values  of  J,  his 
result  involves  a very  large  number  of  parameters  that  are  irrelevant  for 
our  purposes  and  make  it  unsuitable  for  calculation.  In  Appendix  E,  an 
analysis  is  given  in  which  we  sum  over  all  but  one  of  the  unwanted  parameters 
and  obtain 


[pn/J] 

P(J)  = X 
i = l 


(7) 


for  the  probability  of  finding  a crack  of  size  5J  in  a sample  of  n grains  when 
the  probability  of  cracking  for  each  grain  is  p = i - q. 

This  result  may  be  compared  with  a theory  that  neglects  crack  inter- 
action, which,  in  this  case,  would  yield 

Po(J)  = 1 - (1  -pJ)n  (8) 
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A numerical  comparison  of  these  two.  formulas  for  a sample  of  size  n = 100 
and  various  crack  sizes  J is  shown  in  Fig.  2.  It  is  evident  that  the  agreement 
is  very  poor  for  large  or  even  moderate  values  of  p. 

In  Appendix  E,  an  asymptotic  formula  is  derived  for  the  exact  expression 

J -1 

(7)  in  the  limit  n -»  ® , provided  Jp  <?C  1 


HI.  PLANS  FOR  FY  76 


It  is  planned  that  most  of  the  investigations  discussed  in  Section  II  will 
be  completed  before  30  June  1975;  those  that  are  not  will  be  completed  in 
FY  7 6.  It  is  planned  that  the  identification  of  relevant  published  technical 
papers  will  be  completed,  appropriate  data  extracted,  and  the  data  processing 
started  during  FY  75.  However,  it  will  probably  be  many  months  before 
all  the  useful  information  is  extracted  from  the  large  volume  of  data  con- 
tained in  the  world's  technical  literature. 

A.  INCLUSION  OF  POROSITY  AS  A PARAMETER 
IN  THE  THEORY 

This  effort  has  been  started  for  small  pores  in  graphite  (Appendix  C). 
The  approach  followed  is  to  assume  that  a pore  is  a missing  grain  and  that 
then  the  load  is  shared  equally  by  the  remaining  grains  in  the  crack-like 
oriented  array.  For  large  pores,  stress  concentration  factors  must  be 
considered.  In  effect,  the  neighboring  grains  will  be  subjected  to  a greatly 
increased  stress  and  this  will  affect  the  statistics  of  failure. 

B.  REFORMULATION  OF  THE  THEORY  TO  INCLUDE 
THE  STATISTICAL  CONSEQUENCE?  OF  MATERIAL 
VARIABILITY 

This  task  is  of  great  practical  importance.  In  normal  engineering 
practice,  if  the  material  properties  must  be  highly  uniform,  the  rejection 
rate  will  be  high  and  costs  will  increase.  Even  in  the  laboratory,  it  is  nearly 
impossible  to  obtain  identical  specimens,  and  material  variability  obscures 
the  effects  xinder  investigation. 

C.  EXTENSION  OF  THE  THEORY  TO  INCLUDE 
PREDOMINANTLY  COMPRESSIVE  STRESSES 

This  task  is  much  more  difficult  than  the  tensile  stress  task.  In  brittle 
materials  subjected  to  tension,  crack  growth  promotes  instability  and  failure. 
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In  compression,  the  crack  grows  in  such  a direction  as  to  alleviate  the 
stresses,  and  crack  growth  stops.  Thus,  material  failure  in  compression  is 
not  a weakest -link  phenomenon,  but  the  end  result  of  a lot  of  accumulated 
damage.  Because  of  this,  it  is  uncertain  whether  or  not  an  adequate  theory 
can  be  developed  by  the  end  of  FY  76.  However,  the  prospects  are  good  be- 
cause a large  volume  of  new  experimental  information  on  compressive  failure 
of  alumina  is  being  generated  by  M.  Adams  under  the  direction  of  G.  Sines  at 
the  Univ.  of  California,  Los  Angeles. 

D.  EXTENSION  OF  THE  THEORIES  FOR  APPLICATION 

TO  POLYAXIAL  STRESS  CONDITIONS 

For  materials  with  preexisting  cracks,  in  the  approximation  in  which 
it  is  assumed  that  only  the  component  of  stress  normal  to  the  plane  of  the 
crack  contributes  to  fracture,  the  polyaxial  case  has  already  been  treated. 

In  fact,  however,  the  shear  stress  also  contributes;  thus,  the  above  approxima- 
tion is  unconservative.  For  a Griffith  crack  (elliptical  cylinder  stressed  in 
the  plane  normal  to  the  cylinder  axis),  the  contribution  of  the  shear  stress  to 
fracture  can  be  evaluated  exactly  and  the  fracture  statistics  revised  accordingly. 
The  objective  for  FY  76  is  to  go  a step  beyond  this,  if  possible,  and  develop 
a statistical  treatment  for  penny-shaped  cracks. 

E.  REFINEMENT  OF  THEORIES  FOR  GREATER 

ACCURACY 

This  is,  in  a sense,  an  open-ended  task.  Refinements  to  include  porosity 
and  effects  of  shear  have  been  identified  earlier.  It  is  expected  that  additional 
work,  extending  well  beyond  the  work  described  in  Appendix  C,  will  be  done 
on  the  graphite  model  of  Buch,  Zimmer,  and  Meyer.  Variable  grain  size  will 
be  considered.  To  the  extent  possible,  consideration  will  be  given  to  the 
effects  of  crack  interaction  and  variations  in  crack  shape. 
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APPENDIX  A.  SAMPLE  COMPUTER  SEARCH 
REPORT  AND  DIGEST 


A sample  page  of  the  computer  search  of  library  sources  is  shown 
below.  This  is  followed  by  sample  pages  of  the  digest  of  a paper  that  con- 
tains statistical  fracture  data. 


123943a  Fracture  in  polycrystalline  ceramics.  Coble,  R. 
I,.;  Parikh,  N.  M.  (Dep.  Metall.  Slater.  Sci..  Massachusetts 
InM.  Technol.,  Cambridge,  Mass.).  Fracture  1972.  7,  243-314 
(Eng).  Edited  by  I.iebuwitz,  Harold.  Academic:  New  York.  N. 
Y.  A review  is  given  on  the  fracture  in  polycryst.  ceramics,  the 
plastic  and  brittle  behavior  of  ceramic  oxides,  fracture  in 
polycryst.  AljOs,  and  factors  affecting  strength.  130  Refs. 


9602dv  Strength  und  fracture  properties  of  glass-ceramics. 
Hine,  P.;  McMillan.  P.  W.  (Dep.  Phys..  Univ.  Warwick, 
Coventry.  Engl.).  J.  Mat".  Sei.  1973,  8(7).  1041-8  (Eng). 
The  effect  of  an  isothermal  heat  treatment  on  the  strength  and 
fracture  surface  energy  of  a glass-ceramic  derived  from  a Li 
silicate  glass  was  studied.  The  strength  and  effective  surface 
energy  for  crack  initiation  increase  ns  the  mean  free  path  in  the 
interc’ryst.  glass  decreases.  The  strengthening  and  toughening 
mechanisms  are  discussed. 


124020c  Evaluation  of  the  brittleness  of  refractories. 
Pisarenko.  C.  S.:  Gogotsi,  G.  A.  (Inst,  Probl.  Proehn,.  Kiev, 
USSR),  Ogneupory  1974,  (2),  44-7  (Russ).  A Criterion,  *.  for 
the  brittleness  of  refractories  is  proposed:  x “ <>/<»  where  if  is 
the  elastic  and  i , the  limiting  strains,  resp.  A material  is  subject 
to  brittle  fracture  when  x is  equal  or  close  to  unity,  x Decreases 
with  increasing  temp.;  the  decrease  is  pronounced  in  refractories 
of  the  aluminosilicate  type  and  negligible  in  SiC.  T.  Ordentlich 


)37844z  Influence  of  specimen  site  and  mode  of  loading  on 
the  fracture  of  graphite  'Marshall,  P,;  Priddlc,  E.  K.  (Berke- 
ley Nuel.  Lab.,  Cent.  Electr.  Generating  Board,  Berkeley/ 
Gloucestershire,  Engl.).  Carbon  1973,  11(0),  027-31  (Eng). 
Predictions  of  graphite  strength  under  tensile  and  bending  loads 
made  by  using  fracture  mechanics  « ere  compared  with  new  exptl. 
data  and  with  catena,  by  the  method  ol  \V.  Weibull  (1939). 
The  methods  are  comparable  for  similar  specimens,  but  the  Wet- 
bull  method  failed  to  predict  the  obsd.  effects  for  large  specimen- 
site  differences.  The  basis  of  the  Weibull  theory  (a  volumetric 
flaw  distribution)  might  be  inapplicable  to  failure  of  materials, 
resulting  from  surface  defeels.  Examn.  ol  graphite  rods  showed 
a defect-site  distribution  across  the  maternal  that  ltd  to  predic- 
tion of  lower  strengths  lor  the  inner  rod  regions  with  respect  to 
the  outer  regions.  Exptl.  results  followed  the  predicted  trend, 
but  further  examn,  ol  the  fracture-mechanics  techniques  must 
be  undertaken  before  it  can  be  recommended  as  a design  method 
for  graphite  materials. 


9$!02q  Effect  of  air  pressure  on  the  brittle  failure  of 
graphite.  Dubrovskii,  K.  E.;  Kis>el.  V.  V.  (Moscow.  USSR). 
Fit.-Khim.  Afe*h.  Afotrr.  1073,  9(1),  33-5  (Russ).  The 

ultimate  tenule  stress  on  static  bending  (»,)  of  graphite  sealed  off 
from  the  pressure  air  increases  linearly  with  the  pressure  <300 
atm.  A steep  increase  of  a,  was  obsd.  at  higher  pressures.  The 
coui-se  of  the  curve  at  <300  atm.  follows  the  Mirolyubov  law  (I, 
L Col'denblat  and  V.  A.  Kopinov,  1968).  If  air  were  in  contact 
with  the  test  specimen,  the  function  of »,  va.  pressure  was  linear 
over  the  whole  pressure  range  (<500  atm)  and  had  a smaller 
slope  than  that  in  the  previous  case;  the  value  of  a,  at  300  atm 
was  lower  by  —30%.  The  weakening  effect  of  the  air  is  due  to 
the  lowering  of  the  surface  energy  of  graphite  resulting  from  the 
absorption  of  air  by  graphite.  11.  Kucova 


9S!09x  Effects  of  porosity  on  fracture  of  aluminum  oxide. 
Coppola,  J.  A.;  Bradt,  R,  C.  IMater,  Sci.  Dtp  . Pennsylvania 
State  Univ,.  University  Park,  Pa.).  J.  Amrr.  Ceram,  Sac.  1973, 
55(7),  392-3  (Engl.  The  fracture  surface  energy  was  detd,  for  1 
pm  grain  site  AhOs  ceramics  contj,  0.5%  MgO  with  porosity 
<87».  The  porosity  was  located  primarily  at  grain  boundaries 
while  the  fracture  topography  was  predominantly  intergranular, 
The  fracture  mode  was  mixed;  -10%  more  was  trar.sgranulsr. 
The  absence  of  a distinct  decrease  In  the  fracture  surface  energy 
with  decreasing  porosity  suggests  a pore-crack  front  interaction. 


V OL.79,1973 


107553a  Hertzian  fracture  of  a lithium  silicate  glass’ an4 
glass-ceramic.  Nadrau.J.S.;  Rao,  A.S.  (Metall.  Dep..  Unle: 
British  Columbia.  Vancouver.  B.  C.h  J.  Con.  Ceram.  See! 
1972,  41.  63-7  (Eng).  The  fracture  energy  of  Li  silicate  glsu 
prepd  as  a homogeneous  glaxi  is  iets  than  that  of  lie 
phase-srpd.  glass  which  is  less  than  that  of  the  glass-ceramic. 
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IS.  ABSTRACT 

The  investigation  described  herein  was  conducted  to  develop  an  analyti- 
cal model  to  accurately  predict  the  biaxial-fracture  strengths  of  brittle 
materials.  Upon  the  basis  of  work  originated  by  Inglis  and  Griffith,  a new 
biaxial-failure  criterion  has  been  developed.  The  new  criterion  is  a flaw- 
based  theory  concerned  with  the  influence  of  spherical  and  ellipsoidal 
discontinuities  upon  the  fracture  strength  of  brittle  materials.  Also 
considered  in  the  new  theory  are  the  sharp  discontinuities  that  constitute 
a special  limiting  condition  of  the  new  theory  which  coincides  with 
Griffith's  criterion. 

An  experimental  investigation  was  conducted  to  demonstrate  the  validity 
of  the  new  biaxial-failure  criterion.  Oue-hundrcd  and  eighteen  porous  slr- 
conia  test  specimens  were  used  for  the  fracture  tests.  These  tests  were 
conducted  for  uniaxial  tension  and  compression  and  for  two  biaxial-stress 
states  in  the  tension-compression  quadrant. 

A careful  analysis  was  wade  to  develop  test  conditions  that  largely  cir 
cumvented  the  experimental  difficulties  encountered  by  previous  investiga- 
tors. Photoelastic  and  strain-gage  studies  were  conducted  to  verify  the 
experimental  techniques;  the  results  showed  that  homogeneous  stress  fields 
were  developed  in  the  gage  area  of  the  test  specimens* 
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The  experimental  data  were  compared  to  the  predictions  of  the  most 
prominent  failure  criteria.  The  limited  data  available  in  the  literature 
tended  to  fit  the  new  criterion  better  than  the  maximum- .iormal-stress,  the 
Coulomb-Mohr , the  Griffith,  or  the  Weibull  biaxial-failure  criteria. 
Throughout  the  entire  compression-tension  quadrant,  the  fracture  stresses 
obtained  from  this  program  were  in  excellent  agreement  with  the  predictions 
of  the  new  criterion;  no  such  agreement  was  obtained  between  these  data  and 
the  other  failure  criteria. 

. Testing  procedures  had  a pronounced  influence  upon  the  test  results. 
When  the  axial-compressive  load  component  was  applied  before  the  hoop- 
tensile  component,  good  agreement  was  obtained  between  the  data  and  the 
new  criterion,  except  in  the  lower  region  of  the  tension-compression  quad- 
rant. In  this  region,  the  experimentally  determined  strengths  were  higher 
than  those  predicted  by  the  new  criterion.  These  differences  could  be 
explained  by  a mechanism  whereby  favorably  oriented  cracks  were  closed  by 
the  compressive  stresses.  Such  closure  would  permit  stresses  to  be 
.transmitted  across  crack  boundaries,  a condition  inadmissible  in  Inglis* 
original  solution. 

The  excellent  agreement  obtained  between  the  predictions  of  the  new 
criterion  and  the  experimental  data  demonstrates  the  validity  of  the 
analytical  model  and  the  reliability  of  the  experimental  techniques.  Such 
agreement  supports  the  basic  hypotheses  used  to  develop  the  analytical 
■ode.l. 
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Figure  11.  — Test  specimen 


Xw. lined  Avtilrc*  . - ...  • 


1 I limit  ■teuhttd  >1111  etictcl  19  Ut  intuit  ennui  KOI  Unit 
* Itil  ittitiiKiuM  teu.  t't  ceett  -1:1  ui  -51  ut  hi  unit  sum- 
c In  Tift  5-1  let  nucct  m r«»  nil. 


■ TABLE  H 

TEST  RESULTS  FOR  BATCH  FIVE 


1 


I 


f 


SPECIMEN 

TEST 

| DIAMETER 

PRESSURE 

AT 

HOOP 

STRESS 

(mo) 

AXIAL 

LOAD 

(IV) 

AXIAL 

STRESS 

(M!) 

SUITABLE  , 

NO. 

NO. 

MW 

WtoW 

MW 

FAILURE 

<»>|) 

FOR  AHALTSES'  *• 

vhiaxial  noor  tension  retrs 


m 

m 


ioo_ 

“101 


103 


Msa 

3.1530 

3.I3S4 

3.ISSS 

3.ISS1 

3.I5S4 


3.0) 14 
3.6)14 
3 0)14 
30)14 
30)1S 

3.0) 15 


3.1SS3  30)14 


1H0 
I86  0 
3D.0 
30)0 
30)0 
300  0 
303.) 


38S) 

31S0 

m 

3138 

3SI0 

3333 

3430 


AvtUf*  . 


33SS 


NorulueO  A*<  J|t  . 


1 0000 


-M  biaxial,  ran  a 


in 

in 

In 

in 

In 

in 

in 


1)3 

KM 

■HIKM 

Trra 

■IM 

El 

— 'tw 

IM 

Yes 

n< 

rjn 

KM 

1MB 

■m 

— 

■9 

Yes 

Ay«»f« 1 

1 I0US  1 

■1W 

Bi;ni 

..-.A 

.Sil  BIAXIAL  rim  - JTiKOiKO* 


m 

KM 

MMM 

■IW 

— 

KIM 

B9 

MUM 

Yes 

IM 

mm 

EE9 

»W 

■HOI 

TIM 

KM 

mvm 

Yes 

142 

mm 

■>w 

■HW 

■I1M 

■tom 

■>*»■ 

Yes 

I« 

KOI 

■HUM 

■TOM 

HM 

— 

12 00  1 

■HIM 

Yes 

Ai 

TIM 

■HUM 

NotmIizH  AvtoteA  • • .........  ......  | 

■EM 



-j.i  biaxial  rests  - f,ussv*i  Atmco  ntsi  ‘ 


| ME— 

EH 

n> 

BI'HM 

■tom 

H2M 

KM 

■XU 

wkmhmkm 

mm 

BHWl 

KM 

■KM 

KM 

km 

B1SI 

IHKSHIH 

mmntm 

mm 

mma 

im 

Mma 

KM 

■tom 

■!*?!■ 

■kh 

mmsmm 

Kl 

■mm 

1HH 

KM 

KM 

KM 

IBM 

BFMI 

WWW 

KM 

WiWWI 

KTM 

WWW! 

BUM 

BBm 

.1,1  BIAXIAL  lists  - COnSIAHt  inm  l»K‘ 


IX 

mm 

km 

■TOM 

KM 

TIM 

KEM 

■HIM 

Yes 

144 

■M 

KM 

KM 

KM 

TTM 

KM 

■1K« 

Yes 

14) 

m 

BUM 

KM 

■EM 

■TOM 

KM 

KZM 

in 

141 

um 

KM 

■1W 

KM 

— 

km 

■HM 

Yes 

— 

. . 

■rrwi 

, , s,  . „ _ _ . 

NcrsiUteO  A*i»te*. « . 

.....  ... 

KM 

. . 

■TOM 

.... 

AIML  C0»r«(JS*0tf  ?<STI 

IV 

CT 

km 

KM 

... 

m 

1 

■E3M 

Yes 

l* 

KM 

TOW 

KM 

mm 

■TOM 

KM 

Yes 

i» 

El 

KM 

TOW 

mm 

...  ..  . 

■a 

■itviM 

Yet 

10 

El 

■atm 

iBW 



. ..... 

KIM 

tvs 

141  . 

n 

■EM 

■riia 

IWWWW 

■££■ 

■»;•;■ 

Yes 

A 

BPM 

...  . v. 

R.uiut.  A^un*. . 

......  .. 

......... 

......... 

. ..  .1.4 

BIaW 

'SttttutiwulKttf  aiAmAttl  u M ».«»{<  inunul  Kej  mm 
k)«4  tOnuiKtiiM  last.  ouKti  ah  -S  | »t  «4  mm  tun 


* S«  IMt  5-10  let  i»K>n  ea  ixA  Ini  • 


134 


< 

! 


A-IO 


APPENDIX  B.  INTEGRAL  EQUATION  SOLUTION 

FOR  N(<r  ) 

' cr 


In  the  theory  for  volume-distributed  cracks, 
probability  of  failure  in  simple  tension  is  given  as 


* 

the  equation  for  the 


P^(tr)  = 1 - exp 


(B  — 1 ) 


This  can  be  regarded  as  an  integral  equation  relating  P^  and  N,  which  can 
be  solved  as  follows: 


ps  = 1 - Pf  = exp 


■vf  (l  - /^7?) 

J r\ 


£ i-dr 
dcr  cr 
cr 


(B-2) 


in  P s((r)  = 


-vj  (l  - 4 ”•„/«■) 


dN  , 

•z dcr 

d<r  cr 
cr 


(B-3) 


Integration  by  parts  yields 


in  Ps  = 


(T  JT 


■V  (l  - 

o ■'o 


N(cr  ) dcr 
' cr'  cr 


(B-4) 


cr 


*S.  B.  Batdorf  and  J.  G.  Crose,  "A  Statistical  Theory  for  the  Fracture  of 
Brittle  Structures  Subjected  to  Nonuniform  Polyaxial  Stresses,  " J.  Appl. 
Mech.  41  (2),  459  (1974). 
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Since  N(o)  = 0,  the  integrated  term  is  zero  at  both  limits  and  therefore 
vanishes.  Multiplying  both  sides  of  the  remaining  equation  by  <^5- and  then 
differentiating,  we  obtain 


3F  (^<nP,W:)  = -T^ 


N(cr)  = - sfcr  in  Pg  (cr) 
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APPENDIX  C.  SOME  APPROXIMATE  TREATMENTS  OF  FRACTURE 
STATISTICS  FOR  POLYAXIAL  STRESS  STATES 


In  his  original  treatment  of  the  statistics  of  fracture,  Weibull  (Ref.  C-l) 
introduced  a two -parameter  functional  form  for  the  relation  between  simple 
tension  and  probability  of  fracture.  He  also  showed  how  to  compute  the  cor- 
responding statistics  of  failure  for  bending,  torsion,  and  other  stress  states 
that  involve  only  uniaxial  tensile  stresses.  These  applications  have  become 
well  known  and  are  widely  used  (Refs.  C-2  through  C-4). 

In  problems  chat  involve  biaxial  or  triaxial  tensile  stresses,  the  situa- 
tion is  more  complicated.  Weibull  gave,  without  formal  proof,  a procedure 
for  treating  such  problems  and  showed  how  to  apply  it  in  some  simple  cases. 
Some  investigators  have  expressed  doubts  concerning  the  rigor  of  Weibull' s 
treatment  of  polyaxial  stress  states  (Refs.  C-5  and  C-6)  and  there  are  indi- 
cations that  these  doubts  were  later  shared  by  Weibull  himself  (Refs.  C-7 
and  C-8).  Since,  in  addition,  rather  tedious  calculations  are  required  for 
each  polyaxial  stress  state,  in  practical  structures  that  involve  continuously 
varying  stress  states  there  is  a natural  tendency  to  use  approximations.  One 
simple  approximation,  which  constitutes,  in  fact,  the  only  technique  suggested 
for  handling  polyaxial  stress  statistics  in  a recent  treatise  on  fracture 
(Ref.  C-9).  is  to  assume  that 

Ps  <V  V «3>  " PS  ^l'  Ps  <^1  Ps  ^ <C-‘> 

V- 

where  Pg  is  the  probability  of  survival,  and  <r^,  o^,  03  are  t^ie  principal 
stresses.  This  method  was  mentioned  earlier  by  Barnett  et  al.  (Ref.  C-5). 
For  brevity,  it  is  referred  to  here  as  the  Barnett-Freudenthal  (BF) 
approximation . 

One  objective  of  the  present  investigation  is  to  explore  the  limitations 
of  this  approach  and  thus  assist  potential  users  in  understanding  the  general 


C-i 


nature  and  approximate  magnitude  of  the  errors  involved.  Another  objective 
is  to  point  out  an  alternative  approximation  that  is  preferred  in  some  appli- 
cations. For  the  sake  of  simplicity  and.br  evity,  interest  is  focused  primarily 
on  biaxial  tension  applied  to  isotropic  materials  whose  fracture  statistics  in 
simple  tension  can  be  described  with  satisfactory  accuracy  by  Weibull's  two- 
parameter  formulation: 


P£  = 1 - exp  [-k  crm] 


(C-2) 


or 


inPg  = in  (i  - Pf)  = -k  <rm 


(C-3) 


where  P^  is  the  probability  of  fracture.  An  exact  treatment  of  the  Weibull 
theory  solution  is  also  given. 

The  BF  approximation  would  be  strictly  correct  if  all  the  cracks  were 
oriented  with  their  planes  normal  to  any  one  of  the  principal  stresses.  How- 
ever, this  is  rarely  the  case  and  certainly  does  not  occur  in  isotropic  ma- 
terials. Normally,  there  will  be  some  cracks  inclined  to  the  principal  axes 
that  will  be  fractured  by  combined  stresses  even  though  capable  of  surviving 
any  of  the  principal  stresses  applied  individually.  As  a result,  Eq.  (C-l)  is 
usually  unconservative. 

In  the  case  of  equibiaxial  tension,  Eq.  (C-i)  reduces  to 


Ps  {«,  <r)  = [Ps  (<r,  o)]2 


or 


In  Pg  (<r,  cr) 
in  Ps  (<x,  o)  = 2 


(C-4) 


(C-5) 


C-2 


t 


Equation  (C-5)  implies  that,  for  the  small  values  of  of  primary  interest 

in  Pg  (cr,  <r)  in  [1  - Pf  (or,  <r)]  Pf  (<r,  <r) 

in  Pg  (c r,  o)  = in  Li  - (<r,  o)J  “ Pf  (<r,  o)  = 2 (C“' 


i.  e. , the  probability  of  failure  in  equibiaxial  tension  is  twice  that  for  uniaxial 
tension. 

As  noted  earlier,  it  is  clear,  from  general  principles,  that  this  is  an 
underestimate  of  the  probability  of  equibiaxial  fracture.  However,  in  order 
to  determine  the  magnitude  of  the  error  involved,  we  must  turn  to  theory. 

For  this  purpose,  we  select  a recent  reformation  of  weakest  link  theory  in 
which  the  flaws  are  identified  as  cracks  (Ref.  C-10),  and  the  polyaxial  stress 
statistics  are  derived  in  a straightforward  manner  from  the  basic  assump- 
tions. In  Ref.  C-8,  it  was  noted  that  when  the  probability  of  fracture  in  sim- 
ple tension  obeys  Eq.  (C-2),  the  probability  for  equibiaxial  tension  becomes 

Pf  (o-,  <r)  = 1 - exp  ^ (2m  + 1)  k<rm  (C-7) 

For  integral  values  of  rn,  Eq.  (C-7)  becomes 


k (m)  '.  <r  

5)  (m  - 1.5)  — 0.5 


Combining  Eqs.  (C-2)  and  (C-8).  we  obtain 


!n  ps  <T-  mj 

In  Fg  (or,  o)  ” (m  - 0.5)  (m  - 1.5)  — 0.5 


(C-9) 


In  Fig.  C-l,  this  result  is  compared  with  the  previous  one  given  in 
Eq.  (C-5),  and  we  conclude  that  the  BF  approximation  leads  to  the  correct 
result  for  equibiaxial  tension  when  m = i.  However,  for  m = 5 and  iO,  it 
underestimates  small  equibiaxial  failure  probabilities  by  factors  of  approxi- 
mately 2 and  3,  respectively.  Note  that  Eq.  (C-9)  can  be  approximated  sur- 
prisingly well  by  the  very  simple  equation 


in  Pg  (<r,  «r)  ^ 
In  Ps  (<r,  o)  ~ 


2 (m) 


0.45 


(C -10 ) 


In  view  of  the  expected  unconservatism  of  the  BF  approximation;  it  is 
somewhat  surprising  that  Fig.  C-l  implies  that  it  is  correct  for  m -•  1 and 
conservative  for  m<  1.  The  explanation  may  be  that  these  are  abnormal 
cases.  With  the  use  of  Eq.  (C-2),  it  is  readily  shown  that  the  slope  of  the 
P^  (<r)  curve  is  zero  at  <r  = 0 when  m>  i,  in  agreement  with  observation. 

For  m = 1,  the  slope  is  finite;  for  m<  1,  it  is  infinite;  such  behavior  is 
rarely,  if  ever,  encountered  in  nature. 

We  now  consider  the  way  the  size  of  the  error  in  Eq.  (C-l)  varies  with 
stress  ratio,  again  focusing  our  attention  on  the  biaxial  case.  It  was  pointed 
out  earlier  (Ref.  G-8)  that,  for  materials  that  obey  Eq.  (C-2),  the  Weibull 
theory  and  that  of  Batdorf  and  Crose  (Ref.  G-iO)  lead  to  the  same  results.  In 
the  general  case  of  triaxial  tension  characterized  by  principal  stresses  cr^,  <r2 
and  cr\j,  where  <Tj  > o"2  > (r^  > o,  both  theories  imply  that  the  probability  of 
fracture  takes  the  form 


2tt 


P,  = 1 - exp 


-k j d <pj  ^CTjCos  0 sin  0 + <r2  sin  0 sin  0 


'O  o 


2 \m 

+ cr_  cos  0 J sin  8 d 0 


(C-ll) 
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m 


Figure  C-l 


. Ratio  of  in  of  Probability  of  Survival  in 
Equibiaxial  Tension  to  that  in  Uniaxial 
Tension  [equal  to  Pf(cr,  cr)/P£(cr,  o)  when 


where  0 and  <p  are  the  polar  and  azimuthal  angles,  respectively.  As  noted 
by  Shur  (Ref.  C-ll)  the  integral  in  Eq.  (C-ll)  leads  to 


( ) 


m 


d0  = 4ir  cr 


m (ml  )2 

1 (2m  + 1)1 


m 


i»  j>  k 


.sj  tk 

(i‘.  j'.  k'.  )2 


(C  — 12) 


where  i,  j,  k are  any  positive  integers  that  satisfy  the  inequalities 
o < i,  j,  k < m and  i + j + k = 1,  and  where 

s = o^/tTj  (C-13a) 

« 

t = <r3/cr1  (C  -13b) 


For  the  biaxial  case,  <r3  = t = o,  and  the  probability  of  survival  becomes 


where 


Pg  = exp 


4uk 
2m  + 1 


F (m,  s) 


F (m, 


s)  a 


lm‘.  r ST'  !2i)  1 f2to  - ill  1 

ti'-  (m-i)'.]2 


(C-14) 


(C  -15) 


Thus,  in  the  general  biaxial  case. 


in  Ps  <r2)  _ F (m,  s) 

in  Pg  (cr  t,  o)  " F (m,  o) 


F (m,  s) 


where  mis  any  positive  integer,  and  o < s = — < 1. 


(C-16) 
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The  BF  approximation  gives,  for  the  biaxial  case 


in  Pg  (o-j,  c r2)  in  Pg  (c^,  s) 

in  Ps  (<rr  o)  = in  Pg  o)  = 1 + S 


( C -17) 


Curves  that  represent  Eqs.  (C-16)  and  (C -17)  are  shown  in  Fig.  C-2. 

It  is  evident  that,  for  values  of  m of  practical  interest,  the  BF  approximation 
consistently  underestimates  the  fracture  contribution  of  the  second  stress 
and  that  the  error  is  largest,  in  absolute  terms,  at  least,  for  the  equibiaxial 
case,  s = i . 

Another  way  to  compare  theory  with  the  BF  approximation  is  to  take 
the  ratios  of  Eqs.  (C-16)  and  (C-17).  Thus, 


in  (Ps^Th  _ F (m,  s) 
in  (PsW  1 + sm 


(C  - 1 8) 


This  result  is  shown  in  Fig.  C-3.  It  is  evident  that  the  ratio  increases  mono- 
tonically  from  s = o to  s = 1.  Moreover,  in  the  range  2 < m<  10,  the  ratio 
is  reasonably  constant  near  s = 1. 

t 

The  first  of  these  observations  suggests  a simple  method  for  obtaining 
an  upper  bound  to  the  failure  probability,  or  a lower  bound  to  the  survival 
probability  of  a biaxially  stressed  structure;  viz.,  it  implies  that,  in  every 
element  of  the  structure, 


in  (Ps)Th  F (m,  s)  . F (m,  1)  ~ 0.45 

in<Ps>BF"  l + sm  ‘ 2 


(C  - 19) 


Since  Eq.  (C - 19)  holds  for  every  element  individually,  it  must  hold  for  the 
structure  as  a whole. 
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Thus,  we  conclude  that 


-in(P8)BF  <-in(ps)Th  <-m0-45in(Ps)Bj.  (C-20) 

For  small  values  of  P^,  this  implies  that 

(Pf^BF  — ^f^Th<m  ^f^BF  (C-21) 

Because  the  curves  of  Fig.  C-3  are  fairly  flat  near  s = 1,  in  the  case 

of  structures  whose  loading  is  nearly  equibiaxial,  the  second  inequality  in 

Eqs.  (C-20)  and  (C-21)  approaches  equality.  In  the  case  of  circular  disks 

with  uniform  lateral  or  concentric  ring  loading,  for  instance,  the  (small) 

probabilities  of  failure  should,  according  to  the  theories  discussed  here,  be 
0 45 

quite  close  to  m ‘ times  the  result  given  by  the  BF  approximation,  which 
yields  a simple  closed-form  solution. 

As  mentioned  earlier,  the  preceding  analysis  is  based  on  the  assump- 
tion that  the  uniaxial  tensile  fracture  statistics  of  the  material  under  consid- 
eration can  be  described  with  adequate  accuracy  by  Weibull' s two-parameter 
equation.  When  the  material  is  better  represented  by  his  three-parameter 
equation,  the  analysis  becomes  more  complicated  and  is  beyond  th**  scope  of 
the  present  paper.  However,  the  general  observation  that  Eq.  (C-l)  is  un- 
conservative still  applies,  and  it  is  to  be  expected  that  the  size  of  the  error 
might  be  substantial,  particularly  when  the  Weibull  parameter  m is  large. 
Finally,  it  should  be  pointed  out  that,  for  materials  that  obey  Eq.  (C-2),  an 
accurate  Weibull -type  treatment  of  structural  problems  that  involve  biaxial 
or  triaxial  stresses  is  not  difficult  when  a computing  machine  is  used  in  con- 
junction with  finite  elements.  This  is  because  closed-form  solutions  to  the 
integral  in  Eq.  (C-ll)  are  obtainable;  these  solutions  are  given  in  Eqs.  (C-12), 
(C-14),  and  (C-15). 
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APPENDIX  D.  GRAPHITE  MODELING 


A.  GRAIN  STRENGTH 

We  assume  the  body  is  composed  of  (1  - P)n  identical  grains  of  graphite 
and  Pn  pores;  where  n is  the  total  number  of  pore  and  grain  sites,  and  P is 
the  fractional  porosity.  Each  grain  is  assumed  to  have  a critical  stress  S£, 
the  stress  that  will  fracture  the  grain  when  applied  in  its  weakest  direction, 
i.e.,  parallel  to  the  c axis.  When  simple  tension  is  applied  in  some  other 
direction,  it  is  assumed  that  fracture  occurs  when  the  component  of  stress  in 
the  c direction  reaches  Sc.  Since  the  stress  in  the  c direction  is  then  <r  cos  9, 
where  <r  is  the  applied  stress,  and  0 is  the  angle  between  the  stress  axis  and 
the  c axis,  it  follows  that  each  grain  has  a "strength"  given  by 

<rs  = Sjcos2  0 (D-l) 


B.  CRACK  FORMATION 

The  model  of  Buch,  Zimmer,  and  Meyer  * assumes  that  there  are  chance 

aggregations  of  grains,  all  with  c directions  approximately  parallel  to  the  axis 

of  tension,  that  occupy  penny-shaped  regions  with  planes  normal  to  the  tensile 

axis.  Such  arrays  will  be  abnormally  weak  and  will  fracture  prematurely; 

thus  cracking  will  occur.  As  stress  increases,  the  number  and  size  of  cracks 

increase  also.  Consider  first  the  case  of  zero  porosity. 

We  assume  here  that  the  capacity  of  an  array  to  resist  cracking  is  given 

by  the  average  strength  of  its  constituent  grains.  If  the  distribution  of  grain 

strengths  in  an  array  of  Nq  grains  is  the  same  as  in  the  material  as  a whole, 

and  the  maximum  angle  between  c axis  and  stress  axis  is  0 , then  the  mean 

°max 

strength  can  be  found  as  follows . 

J.  D.  Buch,  J.  E.  Zimmer,  and  R.  A.  Meyer,  An  Analytical  Micros tructural 
Model  for  the  Fracture  of  Graphite,  ATR-74(7425)-4,  The  Aerospace  Corp. , 
El  Segundo,  California  (28  June  1974). 
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We  designate  the  number  of  grains  with  an  angle  of  inclination  less 
than  0 as  N(0).  Since  the  orientation  is  random,  the  probability  of  a grain 
having  an  orientation  between  tp  and  tp  + dip  is  the  ratio  of  the  corresponding 
solid  angle,  2ir  sin  ip  dip,  to  the  solid  angle  of  a hemisphere,  i.e, , 

P(0)  dp  - sin  <p  dip  (D-2) 

Thus,  the  probability  that  the  grain  is  inclined  by  less  than  0 is 

J -0 1 P-3) 

I P(p)  dip  = I 'sin  ip  dip  = 1 - cos  0 
*'o  *'o 

Since  all  Nq  grains  have  an  orientation  angle  less  than  0max>  the  number  N(0) 
is  given  by 


N(0)  = N 


i - cos  0 


o 1 - cos  0 


max 


P-4) 


Now,  the  average  strength  of  the  cracks  is,  by  definition 


or 


crs(0)  dN 


P-5) 


Inserting  Eqs.  P-1)  and  P-4)  into  P-5)  and  carrying  out  the  indicated  opera- 
tion, we  obtain 


or 


= S /cos  0 
c 


max 


P-6) 
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If  Pc  designates  the  probability  that  any  arbitrarily  chosen  grain  is 

inclined  by  less  than  0 , then,  from  Eqs.  (D-3)  and  (D-6) 

max 

P = 1 - cos  0 = 1 - S /<r  (D-7) 

c max  c 

If  we  assume  that  all  grains  are  randomly  oriented  and  independent  of  the 
orientation  of  their  neighbors,  then  the  probability  that  any  selected  group 
of  N grains  have  an  orientation  angle  less  than  9max  is  given  by 

PN  = Pc  * - V*>N  <D-8> 

In  particular,  Eq.  (D-8)  gives  the  probability  that  the  N grains  in  a penny- 
shaped array  centered  at  any  designated  site  will  have  an  average  strength 
of  tr  or  less.  Since  we  have  assumed  that  the  strength  of  an  array  is  equal 
to  the  average  strength  of  its  constituent  grains,  it  follows  that  the  proba- 
bility that  the  array  will  fracture  when  the  stress  in  the  material  is  <r  is 
given  by 


PN  = (1  -Sc/<r)N  (D-9) 

We  note  in  passing  that  in  deriving  Eq.  (D-8),  attention  was  focused 
exclusively  on  the  grains  comprising  the  array,  and  nothing  was  said  about 
neighboring  grains.  Thus,  the  proper  interpretation  of  the  equation  is  that 
P^  is  the  probability  that  N or  more  grains  have  an  average  stress  of  <r,  i.e. , 
that  the  crack  formed  contains  N or  more  grains. 

Up  to  this  point,  we  have  assumed  that  all  sites  were  occupied  by  grains 
that  are  identical  except  for  orientation.  We  now  consider  porosity. 
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2 

Following  the  theory  of  Buch,  Zimmer,  and  Meyer,  we  assume  that  the 
array  contains  NP  pores  and  N(1  - P)  grains,  where  P is  the  volume  frac- 
tion of  porosity.  The  average  strength  of  the  array  then  becomes 

<r>  = -7  I NP  • 0 + N (1  - P)  S /cos  0 1 

N L c max  J 

= (l-P)S/cos0  (D - 1 0 ) 

c max  ' ' 

The  probability  that  an  array  will  be  cracked  at  stress  <r  is  equal  to  the  prob- 
ability that  each  of  the  N(1  - P)  grains  will  be  stressed  at  the  limit  of  its 
strength,  i.e., 


PN  * [‘  - <Sc/,r>o]N“'P)  P-“> 

where  subscript  G is  inserted  in  order  to  emphasize  that  S and  or  refer  to 

c 

the  critical  stress  of  an  actual  grain  (not  a site  average)  and  the  mean  stress 
applied  to  an  actual  grain,  respectively.  If  refers  to  critical  stress  aver- 
aged over  all  N sites,  and  <rto  the  average  stress  at  all  N sites  (therefore  to 
the  average  stress  in  the  material),  then 

s'  . (1  - P)  (SC)G  • ^ P-12) 

where  <rT  is  the  elastic  limit,  and 

0-  = crG  (1  - P)  (D-13) 


2 


See  Footnote  i. 
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Thus,  P^.  can  be  written 


PN  = f1  ' (rL/<r)]N(1"P)  P-14) 

Now,  in  order  to  determine  the  number  of  cracks  of  size  N or  larger 
in  a material  that  contains  a total  of  nP  pores  and  (1  - P)n  grains,  we  mul- 
tiply the  probability  of  having  a crack  that  contains  N or  more  sites  centered 
at  a designated  site  P^  by  the  total  number  of  sites.  In  order  to  obtain  the 
number  of  cracks  that  contain  N sites  but  no  more,  we  subtract  the  number 
of  cracks  that  contain  N + 1 or  more  sites.  Thus,  M^n,  the  number  of 
cracks  that  contain  N and  only  N sites  in  a unit  volume  of  material  with  n 
sites  per  unit  volume  is 


MNn  = * (PN 


dP_. 

p/  ) = _£ N 

*N+i'  dN 


(D-15) 


C.  CRACK  STRAIN 

According  to  the  model  of  Buch,  Zimmer,  and  Meyer,  the  inelastic 
strain  is  caused  by  additional  extension  of  the  stressed  material  that  results 
from  the  presence  of  cracks.  In  order  to  determine  the  stress-strain  curve, 
it  is  necessary  to  know  what  strain  will  result  from  a given  number  of  cracks 
per  unit  volume  of  size  N uniformly  distributed  throughout  an  infinite  medi- 
um. This  strain  is  equal  to  the  relative  displacement  of  two  points  initially 
a unit  distance  apart  on  a line  parallel  to  the  stress  axis.  It  can  be  shown 
that  a crack  of  radius  r located  halfway  between  the  points  and  with  its 
plane  normal  to  the  stress  axis  will  result  in  a relative  displacement  given 
by  the  expression 


3 (1  + v)  (2  - v ) £_ 
3 r tt72  E 


(D- 16) 


r 

. 


j 
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Now 


irr2=Ng2  (D-i?) 

where  g is  the  length  of  a grain  edge  (assumed  square),  so  the  strain  will  be 
3/2 

proportional  to  N . Since  a calculation  of  the  strain  averaged  over  all  the 
positions  at  which  one  crack  per  unit  volume  might  be  located  is  not  available, 
we  simply  assume  an  unknown  proportionality  constant  K and  obtain 

e™  = K (<r/E)  g3  N3/2  (D-18) 

K is  an  unknown  constant  that  depends  slightly  on  Poisson's  ratio  v . It  will 
be  determined  experimentally  later . 

D.  DERIVATION  OF  STRESS -STRAIN  RELATION 

We  can  now  compute  the  total  crack  strain,  which  is  simply  the  sum  of 
the  contributions  of  all  the  cracks: 


e 


c 


= iC  MNn 
N=i 


N 

e 


(D-19) 


where  N^  represents  the  number  of  sites  in  the  largest  crack  at  the  time  of 
fracture.  Since,  at  fracture,  there  is  only  one  crack  of  size  N^  in  the  entire 
volume,  and  at  lower  stresses,  there  isn't  even  one,  the  ceiling  on  the  upper 
limit  can  be  removed  with  negligible  error.  Thus, 

CO 

cc  * Z MNn  £ <D-20-> 

N=1 
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Substitution  of  Eqs.  (D-14),  (D-15),  and  (D- 18)  into  Eq.  (D-20)  yields 


= X K (t/E)  g3  £ N3/2  (l 

N=1  ' ' 


v N{  1-P) 


This  equation  is  simplified  slightly  by  noting  that 


n = 1/g  •,  or 


3 . 

n g = 1 


Now 


= exp 


-li 


- N(1  - P)  in  (l 


and 


3n' 


\ N(l-P) 


-i 


= (1  - P)  inll  exp  - N(i  - P)  in 


>*“(»  -v) 


By  using  Eqs.  (D-22)  and  (D-24)  and  converting  the  sum  to  an  integral 
obtain 


= K(<r/E)(i-P  )(l-~)  [ N1* 

' ‘ J0.S 


-1 


^ exp  - N(i  - P)  in(l 
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(D-21) 


(D-22) 


(D-23) 


(D-24) 
, we 


dN 

(D-25) 


The  integral  can  be  evaluated  by  making  the  substitution 


X a N ( 1 - P)  In  ( 1 


The  result  is 


ec  = 


K (o~/E) 


-1 


(1  - P)  in  1 - 


1.5 


P 

Xo 


where 


-i 


Xq  = 0.5  (l-P)in  l--i 


Now 


X 1 ‘ 5 e -X  dx  = X1*5  e"Xdx  - JX°*1' 


(D-26) 


5 i"x  dX  (D-27) 


(D-28) 


5 e'X  dx  (D-29) 


The  first  term  on  the  right-hand  side  is  the  gamma  function  F (2.5)  = 1.33. 
In  order  to  establish  an  upper  limit  to  the  second  term  approximately,  we 
note  that 


D-8 


s * 2'5 
dx  = TT5- 


0.2  (i  - P)  in  I i - -y- 


For  a typical  porosity  P = 0. 2,  this  upper  limit  is  le°<;  than  11%  of  T (2.  5) 
for  any  value  of  stress  greater  than  1.5  cr^  and  less  than  4%  for  any  value 
of  stress  greater  than  2 tr^.  For  the  sake  of  simplicity,  therefore,  we 
neglect  the  second  term  on  the  right-hand  side  of  Eq.  (D-29)  and  approxi- 
mate Eq.  (D-27)  by 


1. 33  K <r/E 


( 1 - P)  in  1 


From  Eq.  (D-31),  it  follows  that  simple  graphites  with  the  same  porosity  and 
Poisson  ratio  will  all  obey  the  relation 


const.  (<r/E) 


In  1 


with  the  same  value  of  the  constant  in  Eq.  (D-32)  applicable  in  each  case. 

In  Fig.  D-l,  a stress-strain  curve  for  POCO  graphite  is  compared  with 
the  results  of  Eq.  (D-32).  At  Southern  Research  Institute  (SRI),  E was  de- 
termined to  be  1.93  x 10^  psi,  and  cr^  appears  to  be  approximately  2.4  ksi. 
Thus,  the  constant  can  be  determined  by  a single  point  on  the  stress -strain 
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TENSILE  STRESS, 


SRI  RUN  No.:  2120-T7 

MATERIAL:  POCO  GRADE  AXF-5Q  GRAPHITE 

SPECIMEN  No.:  4TT-1 

TEMPERATURE:  70°F 

LOADING  DIRECTION:  TRANSVERSE 

STRESS  RATE:  10, 000  psi  / min 

SPECIMEN  GAUGE  SECTION:  0.188  in.  DIAM  x 1.20  in.  LONG 
✓ 

/ 

/ 

/ FRACTURE  (7.7,  7.86) 


E = 1.93  x 10°  psi 


0.2  all 
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2.4  ■ ELASTIC  LIMIT  STRESS 


UNIT  AXIAL  STRAIN,  10  " in. /in. 

Figure  D-i.  Tensile  Stress  vs  Axial  Strain  at  70°F  for  POCO 
Grade  AX-5Q  Graphite  Specimen  4TT-1  from 
Billet  4(E2A-30A) 
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curve.  The  point  <r  = 7.8  ksi,  ec  = 3.6  x 10  was  chosen  for  this  purpose, 
and  a value  of  0.  2 was  obtained  for  the  unknown  constant.  The  equation  thus 
becomes 


e 


c 


fn 


0.2  (oVE) 


For  convenience,  the  function 


0.2 


(D-33) 


P-34) 


is  plotted  in  Fig.  D-2.  This  function  is  simply  the  ratio  of  the  crack  strain  to 
the  elastic  strain. 

In  Eq.  (D-33),  it  is  implied  that  the  crack  strain  is  determined  by  ap- 
plied stress  and  two  material  parameters,  Young's  Modulus  and  the  elastic 
limit.  The  latter  is  sometimes  difficult  to  determine  accurately.  There- 
fore, for  purposes  of  checking  the  accuracy  of  the  equation,  it  is  sometimes 
preferable  to  determine  what  value  of  <r^  provides  good  agreement  with  the 
stress -strain  curve  near  the  fracture  point,  and  then  determine  if  this  value 
falls  in  the  range  of  values  that  would  be  estimated  by  different  observers  as 
the  first  appearance  of  nonlinearity  in  the  curve.  Solving  Eq.  (D-33)  for 
or.  , we  obtain 


where  e denotes  elastic  strain, 
e 


(D-35) 
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Figure  D-2.  Crack  Strain  vs  Elastic  Strain  for  POCO 
and  ATJS  Graphite 
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COMPARISON  OF  STRESS -STRAIN  PREDICTION 
WITH  EXPERIMENT 


The  constant  0.2  in  Eq.  (D-33)  was  chosen  to  give  good  agreement 
between  theory  and  experiment  at  the  fracture  stress  for  a particular  POCO 
graphite.  The  calculated  points  designated  in  Fig.  D-l,  however,  indicate 
that  the  entire  stress -strain  curve  is  fitted  very  well  by  the  equation.  Fig- 
ure D-3  is  presented  in  order  to  check  the  accuracy  of  the  equation  against 
another  POCO  graphite  with  different  values  of  E and  c . tr^  was  determined 
by  using  Eq.  (D-35).  It  is  apparent  that  the  agreement  between  theory  and 
experiment  is  again  very  good. 

POCO  is  a relatively  simple  graphite  in  that  it  is  essentially  isotropic 
and  has  approximately  uniform  grain  size  and  pores  approximately  the  same 
size  as  the  grains.  In  Figs.  D-4  and  D-5,  Eq.  (D-33)  is  compared  with  test 
data  for  ATJS  graphite,  which  is  anisotropic  and  has  variable  pore  and  grain 
size.  These  complications  are  not  taken  into  account  in  the  theory  underly- 
ing Eq.  (D-33)  and  are  beyond  the  scope  of  this  report.  The  general  shape 
of  the  curve  appears  to  be  predicted  well,  but,  in  some  cases,  the  elastic 
limit  is  noticeably  different  from  the  value  of  that  gives  the  observed 
crack  strain  at  fracture.  It  is  of  some  interest,  therefore,  that  if  the  con- 
stant in  Eq.  (D-33)  is  changed  from  0.2  to  0.13,  very  good  agreement  is 
obtained  with  experimental  data  on  ATJS.  It  is  possible  that  the  differences 
between  simple  and  complex  graphites  mainly  affect  the  constant  rather  than 
the  basic  form  of  the  equation  that  gives  crack  strain  as  a function  of  stress. 
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* ELASTIC  LIMIT  STRESS 

* THEORETICAL  POINTS 

FOR  <7^  *2.8 


SRI  RUN  No.:  2120-T8 

MATERIAL:  POCO  GRADE  AXF-5Q  GRAPHITE 

SPECIMEN  No..-  4TT-4 

TEMPERATURE:  70°F 

LOADING  DIPv'-V»:  TRANSVERSE 

STRESS  RATti  .,0,000  psi/min 

SPECIMEN  GAIJGi.  SECTION: 

0.18b  ‘n.  DIAM  x 1.20  in.  LONG 


Figure  D-3.  Tensile  Stress  vs  Axial  Strain  at  70 °F  for  POCO 
Grade  AXF-5Q  Graphite  Specimen  4TT-4  from 
Billet  4 (E2A-30A) 
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SRI  RUN  No.:  2139-T14-L16-7 

MATERIAL:  ATJS  GRAPHITE 

SPECIMEN  No.:  S-2-4 

TEMPERATURE:  70°F 

LOADING  DIRECTION;  WITH  GRAIN 

STRESS/ STRAIN  RATE:  10,000  psi /min 

SPECIMEN  GAUGE  SECTION:  0.250  in.  DIAM  x 1.2  In.  LONG 


AXIAL  STRAIN,  103  in. /in. 


Figure  D-4.  Tensile  Stress  vs  Axial  Strain  for  ATJS 
Graphite  Specimen  S-2-4 
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TENSILE  STRESS,  103  psi 


SRI  RUN  No. : T— 23 
MATERIAL:  ATJS 
SPECIMEN  No.:  A20W 

TEMPERATURE:  70°F 

LOAD  DIRECTION:  WITH  GRAIN 
STRESS  RATE:  10, 000  psi  / min 


APPENDIX  E.  THE  THEORY  OF  RUNS  APPLIED 
TO  CRACK  SIZE 

(R.  H.  Huddle  stone  and  D.  C.  Pridmore-  Brown) 


The  theory  of  the  distribution  of  cracked  grains  in  a one -dimensional 
model  that  consists  of  a single  row  of  grains,  each  with  a probability  p of 
being  cracked,  is  identical  with  the  theory  of  runs  treated  extensively  in  the 
statistical  literature.  First,  a brief  outline  is  given  of  those  portions  of  the 
theory  that  are  relevant  to  our  problem,  following  the  approach  and  notation 
of  Fisz.  * Then,  these  results  are  used  to  derive  a formula  for  the  proba- 
bility of  occurrence  of  a crack  of  a given  size  or  greater.  This  formula, 
which  we  believe  to  be  new,  involves  a final  summation  over  a number  of 
terms  equal  to  the  number  of  cracked  grains  divided  by  the  crack  length. 
Although  it  has  not  been  possible  to  perform  this  summation  analytically,  an 
asymptotic  expression  is  derived  for  it,  which  is  then  compared  with  a com- 
monly used  heuristic  formula.  Finally,  some  numerical  comparisons  be- 
tween the  exact  and  heuristic  results  are  presented. 

Consider  a one -dimensional  row  of  n grains,  of  which  n^  = pn  are 
cracked  and  n-,  = n - n^  are  uncracked.  Then,  there  are  ( n^)  ways  in  which 
the  cracked  grains  can  be  distributed  among  the  uncracked  ones.  Here,  a 
run  of  i contiguous  cracked  grains  is  referred  to  as  a crack  of  length  i,  and 
the  total  number  of  such  cracks  of  length  i is  denoted  by  kj..  Similarly,  the 
number  of  runs  of  uncracked  grains  of  length  i is  denoted  by  k£..  Finally, 
the  total  number  of  cracks  of  all  lengths  is  denoted  by  kj,  and  of  noncracks 
by  k£-  Then  ‘ . 


Skli  = kl  ' SkEi  = k 

Eikj.  = n^  ; ^*k2i  = n 


2 

2 


(E-i) 


Fisz,  Probability  Theory  and  Mathematical  Statistics,  John  Wiley, 
N ew  Y ork  (1963). 


E-i 


where  the  summations  are  from  1 to  a on  the  left  and  1 to  n 2 
on  the  right. 


Of  the  total  number  of  configurations  or  ways  of  distributing  n^  cracked 
grains  among  n2  = n - n^  uncracked  ones,  vix. , that  number  that  con- 

tains exactly  k^  cracks  of  length  1,  kj2  of  length  2,  . . . k^n  of  length  n^, 

k0  _ noncracks  of  length  n0  . . . is  just 


k9.  noncracks  of  length  1, 


noncracks  of  length  n- 
2n2  6 2 


n(i 


ii'  *2i» 
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where 


denotes  the  multinomial  coefficient 


kj 


V ' k 1 
KU  * Ki2  * 


"In. 


which  represents  the  number  of  ways  that  k^  objects  can  be  distributed  among 
n^  piles  such  that  k^  objects  are  in  the  ith  pile,  and  G(x)  = 2 for  x = 0,  1 for 
x = 1,  and  0 otherwise.  G ^ |k^  - k2  accounts  for  the  fact  that  the  cracked 
and  uncracked  runs  must  alternate  such  that  kj  and  k2  can  differ  at  most  by 


one. 


On  the  assumption  that  each  of  the  ( arrangements  of  cracked  and 
uncracked  grains  is  equally  probable,  the  ratio  of  Eq.  (E-2)  to  (n^becomes 
the  conditional  probability  that  tlvre  are  kj.  cracks  of  length  i,  i = 1,  n^, 
and  k2j  noncracks  of  length  j,  j = 1,  n2,  given  that  ther«  are  n^  cracked 
grains  and  n2  uncracked  grains 


P (k 


li’  ~2i  ' ni 


n. 


(E-3) 
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Now,  we  compute  the  probability  that  there  is  at  least-one  crack  of  a given 
size,  say  J,  or  greater.  Thus,  we  take  Eq.  (E-3),  which  contains  n param- 
eters (n^  kjs  and  ^k^s),  and  sum  over  all  the  k^s  and  then  sum  sepa^^cely 
over  all  the  k^'s  for  i < J and  all  the  k^  for  i > J such  that  we  end  up  with  an 
expression  for  the  probability  that  contains  the  single  parameter  J.  The  first 
of  these  steps  was  carried  out  by  Fisz  with  the  identities 


- 1 
- 1 


and 


which  are  straightforward  to  derive.  Thus, 


N 


At  this  point,  Fisz  summed  over  the  k^  to  obtain 


N(; 


:i*  nl’ 


(E-4) 
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which  is  the  number  of  configurations  with  exactly  cracks  of  any  length. 
Instead,  we  sum  over  the  long  and  the  short  cracks  separately.  Put 


i < J 
i > J 


with 


2 s 


li 


= s 


1' 


s ia  - f1# 


!1+*1 


k. 
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Then,  it  is  clear  that 


Thus, 


(E-5) 


We  now  sum  Eq.  (E-5)  over  s^.  andi^..  Note  first  that 


+ x2  + 


+ x 


J-i 


can  be  expanded  on  the  one  hand  in  a multinomial  expansion 


S1 

sii 


2is 

x 


ii 
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where  the  summation  over  all  sets  is  s^.  such  that  2s^  = s^,  and  on  the 
other  hand  in  a binomial  expansion 


•i/l-x'-V  SI 
x i — = x 


1 /sJ 


(‘-x)'s‘E(.7 (- 1,1 


(j-i)i 


i=0 


Similarly 


li 


Sii , . Ji 
x 11  = x 


thus 


(x  + *2  + . . .+xJ-1)S‘(xJ+xJ'+1  + . . .)'> 


■E 


sa 


Sis,. 


'ii 


2uii 


S,+J1, 


= X 


si+Jii 


(J-l)i 


m=0 


)xm-i:(i1)<-i>ix<j-i,i 

i=0 
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I 


We  now  require 


Sis  . . + Sii, . = Sik. . 
li  li  li 


to  equal  in  the  preceding  expression.  Then 


sh 


si  *i 


when  the  summation  is  over  all  partitions  |s^|i  = 1,  . . . J-l|  of  ^ n j* 

|i  = J,  . . . n^  of  is  equal  to  the  coefficient  of  x in  the  preceding 
expression.  This  coefficient  is  composed  of  those  terms  in  the  summation 


for  which 


s1  + Jf1  + m+(J-  l)i  = n^ 


nl  " S1  " Jil  - (J  " A)  1 


Thus 


nt  - 1 (J  - 1)  (1,  + 1 
B1 + 11  • 1 
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or,  equivalently 


kj,  Sj,  nj,  n2 


)“ 


kl  ' 1 


(-  l)1 


It  does  not  appear  possible  to  perform  the  summation  over  i.  However, 
a summation  over  kj  can  be  carried  out  as  follows.  First,  repeated  use  of 
the  binomial  identity 


- C 

- C 


permits  one  to  reduce  the  number  of  terms  that  contain  kj.  Thus 


n. 


+ 


i - 


n/lj+l-i-Jj 
/U,  -<i  -i 
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Thus 


N 


(kr  iv  nv  n2)  = ^<-  l)1^. 


n,  + 1 \ / n,,  + 1 - i \ 


2 ‘ * \ / “2 
*1 


kl  ~*l  ~l 


nj  - i - (J  - i) 
ki  ‘ 1 


(i1  +i)j 
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We  now  sum  over  by  using  the  identity 


Note  that  can  run  from  kj  = 1 to  kj  = n2  + 1,  which  is  the  upper  limit 
obtained  when  each  uncracked  grain  acts  as  a crack  separator.  Thus 


N nr 


n. 


n2+1_il 


)■  E (-«* 


i=0 


n - J (I1 
n2 
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If  we  now  sum  over  all  f starting  from  i ^ 
upper  limit,  we  expect  to  find 


= 0 and  going  to  the  natural 
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Thus,  it  is  concluded  that 


I 


d 


n0+i\  /n2+i-i\/n-J(ij+i) 


(n_  f i\  t n„ 

• E - 

i=u  ' je±=o  x 1 
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and  hence 


N ( > 1,  nj,  n2  j = N | n^  n2  ) - N ( nz) 
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n2  + 1 ^ / n - Ji , 
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i=0 
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This  expression,  which  is  henceforth  denoted  by  N (J,  n^,  n2),  represents 
the  number  of  configurations  that  have  at  least  one  crack  of  length  J or  greater. 


f ■ 

$ 

* 


We  now  assert  that 


n2+1 


,n2  + 1\  , n - Ji^ 


Ehi‘1 


n. 


= 0 


i=0 


This  can  be  proved  by  using  the  identities 


I i 


i- 


E-9 


and 


n 


E<-‘»n+i  (. 


t n \ . / p + qi ' 


n 


= q 


i=0 


the  second  of  which  has  been  proved  by  Riordon.  From  this  r.esult,  it  fol- 
lows that  expression  (E-10)  is  equal  to  ( ) when  J < [n/n^^^]  and  that, 

otherwise,  the  upper  limit  on  the  summation  may  be  taken  to  be  [n^/J], 
where  the  square  brackets  denote  integer  part  of. 

Dividing  through  by  ( ),  we  find,  for  the  probability  of  at  least  one 

M 

long  crack  (of  size  > J) 


= 1 for  J < [nj/n£  + 1] 


(E-ll) 


The  corresponding  unconditional  probability  is 


P 


n2 


) = P(j|nr 


n 

(1  -P) 
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Strictly  speaking,  this  expression  should  be  summed  over  n^.  Instead,  we 
believe  it  is  sufficient  to  retain  expression  (E-il)  and  merely  replace  n^  and 
n2  with  their  expected  values  pn  and  qn,  respectively,  where  q = 1 - p. 


2 


J.  Riordan,  Combinatorial  Identities,  John  Wiley,  New  York  ( 1968),  p.  51. 
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We  now  derive  a closed-form  analytical  expression  for  Eq.  (E-il) 
that  is  valid  asymptotically  as  n -*  i.e. , for  large  samples. 

First,  expression  (E-ll)  is  rewritten  as 

p <J)  (nq  + l)li)  qno  (‘  "rb ) 

which  follows  directly  from  the  definition  of  the  binomial  coefficients.  Here 


(nq+  l)(i)  = 


2-i 

n (nci + k) 

k=l 


= exp  2 in  (nq  + k) 


= (nq)1  [1+0  ( 1/n)] 
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Simularly, 


nq  - i 

fl  (‘ -^?k)  =exPEin  (‘ -rh) 

k=0 


Thus,  as  n-*“ 


P(J)- 


[pn/J] 
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Expression  (E-13)  represents  a finite  sum  of  terms  of  alternating  sign. 
For  large  n,  the  partial  sums  oscillate  with  increasing  amplitude  until  a 
value  of  i is  reached  that  is  comparable  to  enc*  ^ . After  this,  the  ratio  of 
successive  terms  is  less  than  unity,  and  the  partial  sums  converge  rapidly. 
Thus,  if  the  upper  limit  [pn/J]  is  much  greater  than  this  critical  value  of  i, 
i . e . , if 


pn/  J » e11^  ^ 

(E-14) 

qJ  PJ-1  « i/e 
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LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  Is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the'many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertised  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory':  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared. technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Physics  Laboratory:  Atmospheric  and  ionospheric  physics,  radla. 
tlon  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  alfglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  In  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  apace  cn  space  systems. 


THE  AEROSPACE  CORPORATION 
El  Segundo,  California 


